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We present a new action for highly improved staggered fermions. We show that perturbative calculations for the new action 
are well-behaved where those of the conventional staggered action are badly behaved. We discuss the effects of the new terms 
in controlling flavor mixing, and discuss the design of operators for the action. 



1. INTRODUCTION 

Naive and staggered fermion actions (which 
have very similar properties) have numerous advan- 
tages over actions based on four-component Wilson 
fermions. They have an exact chiral symmetry, al- 
lowing a simpler approach to the light-quark limit. 
They require only even dimension correction oper- 
ators, simplifying the improvement program. Stag- 
gered fermion propagators can be calculated much 
more rapidly than Wilson fermion propagators. 

They thus offer the possibility of performing much 
higher accuracy calculations on current computers 
than with other methods. However, they have his- 
torically suffered some drawbacks which must be ad- 
dressed before they can be used effectively. Glu- 
ons of momentum K induce unphysical transitions 
between different flavors of doubler quarks, leading 
to large flavor symmetry breaking, especially in the 
pion sector. The perturbation series, even using im- 
proved perturbation theory, is often more poorly be- 
haved than the series for Wilson fermions. With stag- 
gered fermions, quarks of different spins live on dif- 
ferent sites of the lattice, leading to large O (a^) errors 
that are not present for Wilson fermions. Finally, the 
doubling of the fermion spectrum requires the taking 
of a root of the determinant in simulations, requiring 
more care in understanding the locality of the theory. 

This paper will show that the improvement pro- 
gram provides an elegant solution to the first three of 
these problems. 



2. IMPROVED STAGGERED FERMIONS 

One significant problem with ordinary naive and 
staggered fermions is the large flavor nondegeneracy. 
This problem is worst in the pion sector. For the B me- 
son, and other particles with a single light quark, there 
is no flavor mixing at the valence level, because the 
heavy quark would be driven far off-shell by a gluon 
of momentum K. However, since Bs couple to virtual 
pions, flavor mixing interactions can still be impor- 
tant. 

At tree level, they arise from transitions between 
doubler quarks of different flavors induced by gluons 
of momentum K In Ref. [Bh, it was shown how 
to turn "fat-link" improvement[Q|, which suppresses 
the coupling of quarks to these gluons, into a tree level 
O (fl^) improved action by proper tuning of the coeffi- 
cients and the inclusion of an additional term. 

The improved action for naive quarks is given by: 

5imp = Y,w{x)(^y-A'-jy-A^ + m^\v{x). ( 1 ) 

(The connection between naive and staggered quarks 
is the usual one [|^]). is a smeared link variable, 
identical to the ordinary link variable (up to eiTors of 
(fl^)) for low gluon momentum, but vanishes when 
a gluon of momentum n/a is extracted. The action 
has an exact doubling symmetry 
\\f(x) B^(x)\(/(x), ^{x) where 

^t^i^) = np ('YsYp) exp(/^7t -x/fl) for one of fifteen 
vectors ^ with one or more non-zero components. 
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C= (1,0,0,0), (1,1,0,0), etc. 

3. HIGHLY IMPROVED STAGGERED 
FERMIONS 

Monte Carlo calculations showed a large reduc- 
tion in pion flavor nondegeneracy with this action [^]. 
However, significant flavor breaking still remains, so 
further improvement is desirable for high precision 
calculations. In Ref. [|[) a nonperturbative determi- 
nation of gluonic corrections to the tadpole improved 
tree level o(fl^) improved staggered action was per- 
formed. It found a small improvement in pion fla- 
vor breaking was possible, but not the large reduction 
that is still desirable. Therefore, incorporation of ad- 
ditional operators into the action is needed. 

The operators required for full o(ajfl^) improve- 
ment are four fermi contact interactions [^. They 
arise from transitions in quark-quark scattering be- 
tween quarks of different doubler flavor induced 
by two gluon exchange. These contact terms can 
be transformed into interactions between fermion- 
bilinears and auxiliary scalar fields. The possible op- 
erators are governed by chiral and doubling symme- 
tries. 

The required Lagrangian for the quark-scalar field 
interactions is 



Jt{x)^%{x)+h,{x)^s,{x), (2) 
, 8 is a color index, and where the <t> 



where b = \ 

fields are expressed in terms of normalized auxiliary 
scalar fields (|) according to 

etc. The projection operator T (Q is given by 



{^,)'P{Q^>l^{x\ 



(3) 



2'(Q=n 



(4) 



An independent Gaussian random distribution is gen- 
erated for each component of each auxiliary field. 

The color-octet vector current f^{x) for staggered 
fields is given by 

^ V=±)l 



X{x)t'u^{x)x{x- 



flV 



^X{x + av)U^ix)T'xix) 



(5) 



where r|^(ji:) is the usual staggered phase. 

The color-octet axial-vector current for stag- 
gered fields is given by 



JU^) - 17^A<WL xix)T'll{x;x + p)xix + p) 
+ X{x + p)uHx;x + p)T''xix) 



(6) 



where the sum over p is over the eight vectors of 
length \/3a perpendicular to fi 



p2 = 3fl^ pjj^O, 



(7) 



and 

U (x;x- 



P) = 



Upj{x)Up2{x + api) 
xUp^{x + api+ap2] 



symm. 



(8) 



where the symmetrization is over the six assignments 
of the components of p to the indices pi, p2, and 
P3. The symmetrization is done so as to avoid in- 
troducing new discretization errors. The phase T^(x) 
comes from diagonalizing the axial vector current, 
\\l{x)j5jf,\\i{x + p). 

Finally each of the coefficient functions, is ex- 
pressed in terms of two scalar functions of 



c,c|f'(c^)+(i-c,)cr(c'), 



(9) 



and analogously for Cy ' , c^'f^ , and c''^^ . The results are 



given in Table |l]. 



4. OPERATOR DESIGN: NAIVE VS. STAG- 
GERED FERMIONS 

Because of the transformation that turns naive 
fermions into four identical copies of staggered 
fermions, it is clear that a simulation with naive 
fermions (taking the sixteenth root of the determi- 
nant) and one with staggered fermions (taking the 
fourth root) would give identical results. Staggered 
fermions are therefore clearly preferable in simula- 
tions since they are a factor of four cheaper. For va- 
lence quarks and their operators, on the other hand, 
we should examine each for advantages and disad- 
vantages. For staggered fermions, different spins live 
on different lattice sites. Such operators as the (a^) 
improved four-quark operator used in Bk are spread 
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Table 1 

One-loop coefficients for the fermion bilinears in 
Eq. (^), using the Symanzik-improved gluon action 
to compute the gluon propagator. The results are in 
units of ay (^* w Ti/a), and are accurate to a few parts 
in the last digit. Entries with an / are imaginary. 





''V 








1 


0.880 i 


0.500 


0.643 i 





2 


0.435 i 


0.438 


0.217 i 





3 


0.335 i 


0.409 


0.244 i 





4 


0.300 i 





0.220 i 







•-A 


-m 

'-A 


•-A 




1 


0.404 i 


0.518 


0.295 i 





2 


0.364 


0.228 i 





0.166 i 


3 


0.198 i 


0.198 


0.145 i 





4 


0.190 












over 4^ blocks on the lattice. It is reasonable to ex- 
pect large {ap)^ discretization errors from such oper- 
ators, as are indeed observed. For naive quark opera- 
tors such as (\(7(x)ri(/(x))^, all quarks live on the same 
site. Hence, at tree-level they have no finite a errors 
to all orders in a. 



5. CONVERGENCE 
THEORY 



OF PERTURBATION 



Table 2 

For improved staggered quarks (here with Wilson 
glue), unlike unimproved staggered quarks, all coef- 
ficients in one-loop operator renormaUzations are of 
order one. 





Unimproved 


Improved 


Operator 


Staggered 


Staggered 


\j/V|/ and Yii^l 


l-4.17ai/ 


l-0.39av 


iij>Y"\|/ and Ytl^'^ 


l-1.57av 


l+0.48av 




1 -0.56av/ 


1 +0.66av 


VjAj/0+l|A|/ 


l-1.76av 


l + 1.26av 


\jAj/0_l|A|/ 


1 -5.57av 


l+0.40av 


ijAj/Oll|A(/ 


1 -2.40av 


l + 1.12av 


\jAj/02W 


1 -8.11av 


l-0.17ai/ 



theory program of Ref. [|8|] does not work well is 
staggered fermions. This was later explained by the 
existence of fermionic tadpole diagrams in naive and 
staggered fermions: |Jl|,^ fermion loops from all dou- 
bler quark flavors, created by gluon loops of momen- 
tum 71. Since improved staggered fermions suppress 
precisely the coupling to these gluons, we should ex- 
pect that these effects are absent here and that per- 
turbation theory is well behaved. Table ^ contains 
one-loop renormalization coefficients for some of the 
naive fermion operators advocated in the previous 
section. It is clear that for the improved case, the co- 
efficients of a are all of order one, as argued, while in 
the unimproved case they range as high as eight. 

6. CONCLUSIONS 

We have shown that the improvement program ele- 
gantly removes many of the traditional difficulties of 
naive and staggered fermions. In future work, we will 
expand on the work presented here, and apply the new 
action in Monte Carlo calculations. 
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The one case in which the improved perturbation 



